Generating spin-entangled electron pairs in normal conductors using voltage pulses 



in 
o 
o 

(N 
•*-> 

o 
O 

00 



A.V. Lebedev a , G.B. Lesovik Q , and G. Blatter 6 
1 L.D. Landau Institute for Theoretical Physics, RAS, 119334 Moscow, Russia and 



Theoretische Physik, ETH-H6nggerberg, CH-8093 Zurich, Switzerland 
(Dated: February 2, 2008) 

We suggest an operating scheme for the deliberate generation of spin-entangled electron pairs in 
a normal-metal mesoscopic structure with fork geometry. Voltage pulses with associated Faraday 
flux equal to one flux unit <E>o = hc/e drive individual singlet-pairs of electrons towards the beam 
splitter. The spin-entangled pair is created through a post-selection in the two branches of the fork. 
We analyze the appearance of entanglement in a Bell inequality test formulated in terms of the 
number of transmitted electrons with a given spin polarization. 
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I. INTRODUCTION 

Quantum entanglement of electronic degrees of free- 
dom in mesoscopic devices has attracted a lot of inter- 
est recently. Early proposals for structures generating 
streams of entangled particles exploit the interaction be- 
tween electrons as a resource for producing entanglement, 
the pairing interaction in superconductors^ or the repul- 
sive Coulomb interaction in confined geometries^. Re- 
cently, another class of devices has been suggested which 
avoids direct interparticle interaction; instead, the en- 
tanglement originates from a proper post-selection of or- 
bital^££ or spini&S degrees of freedom. The majority of 



these proposals deals with the situation where the entan- 
gled particles are emitted in a random and uncontrolled 
fashion, while entanglement on demand is implicit in the 
scheme proposed in Ref. 

Current interest concentrates on setups which are ca- 
pable to produce pairs of entangled electrons 'on de- 
mand'. Such controlled entanglement is an essential step 
towards the realization of quantum computing devices 
for which electronic orbital- or spin degrees of freedom 
may serve as qubitsiS. In addition, prospects to convert 
electronic entanglement into a photonic one with high 
efficiency look promisingii; this may open new opportu- 
nities for the manipulation of entangled photons with an 
enhanced efficiency. 

While entanglement on demand is implicit in the work 
of Ionicioiu et ai, a detailed discussion of the controlled 
production of entanglement in a mesoscopic device has 
only been given recently by Samuelsson and Biittikcr 12 ; 
they proposed a scheme for the dynamical generation 
of orbitally entangled electron-hole pairs where a time- 
dependent harmonic electric potential is applied between 
two spatially separated regions of a Mach-Zehnder inter- 
ferometer operating in the Quantum Hall regime. This 
(perturbative) analysis concentrated on the limit of a 
weak pumping potential generating only a small fraction 
of entangled electron-hole pairs per cycle. Later on, sev- 
eral schemes have been suggested producing entangled 
pairs on demand with a high efficiency: in their setup, 
Beenakker et al. 13 make use of a ballistic two-channel 
conductor driven with a strong oscillating potential. In 



their non-perturbative analysis they demonstrate that 
this device can pump up to one (spin- or orbital) entan- 
gled Bell-pair per two cycles. A different proposal based 
on spin resonance techniques acting on electrons trapped 
in a double quantum dot structure and subsequently re- 
leased into two quantum channels has been suggested 
by Blaauboer and DiVincenzoii; their detailed analysis 
of the manipulation and measurement schemes demon- 
strates that the production and measurement of entan- 
gled pairs via an optimal entanglement witness can be 
performed with present days experimental technology. 



up lead 




FIG. 1: Mesoscopic normal-metal structures in a fork geome- 
try generating spin-correlated electrons in the two arms of the 
fork. The Bell type setups detect the number of transmitted 
particles Ni, i = 1,3, with spin projected onto the directions 
±a in the upper arm and correlates them with the number Nj , 
j — 2, 4, of particles with a spin projected onto the directions 
±b in the lower arm. (a) Fork with a simple splitter with 
particles injected along the single source lead, (b) Fork in the 
geometry of a four-terminal beam splitter with particles in- 
jected along one of the incoming channels only. Quenching the 
transmission T u< j between the upper and lower leads allows to 
eliminate equilibrium fluctuations spoiling the entanglement. 
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In the present paper, we discuss an alternative scheme 
generating pulsed spin-entangled electron pairs in a 
normal-metal mesoscopic structure arranged in a fork ge- 
ometry, see Fig.^ In this device, spin-entangled electron 
pairs are generated via the injection of spin-singlet pairs 
into the source lead from the reservoiA This entangle- 
ment is made accessible by splitting the pair into the 
two leads 'u' and 'd' and subsequent projection (through 
the Bell measurement) to that part of the wave func- 
tion describing separated electrons travelling in different 
leads&£. Rather the n quan tum pumping with a cyclic 
potential as in Refs. Il'2ll.l our proposal makes use of 
definite voltage pulses generating spin-entangled electron 
pairs. A pulsed sequence of ballistic electrons is im- 
plicitly assumed in the generation of orbital entangle- 
ment by Ionicioiu et al£, however, no description has 
been given how such (single-electron) pulses are gener- 
ated in practice. Below we discuss a scheme where volt- 
age pulses of specific form accumulating one unit of flux 
$o = —c J dt V(t) and applied to the source lead V gen- 
erate pairs of spin-entangled electrons which then are dis- 
tributed between the two outgoing leads of the fork, the 
upper and lower arms denoted as 'u' and 'd'. These spin- 
entangled electron states are subsequently analyzed in a 
Bell experiment 15 involving the measurement of cross- 
correlations 16 between the number of electrons transmit- 
ted through the corresponding spin filters in the two arms 
of the fork, see Fig. ^ Using time resolved correlators, 
we are in a position to analyze arbitrary forms of voltage 
pulses and determine the resulting degree of violation in 
the Bell setup. We find that Lorentzian shaped pulses 
generate spin-entangled pairs with 50 % probability, cor- 
responding in efficiency to the optimal performance of 
one entangled pair per two cycles as found by Beenakker 
et ali&. The reduction in efficiency to 50 % is due to the 
competing processes where the spin-entangled pair gen- 
erated by the voltage pulse propagates into only one of 
the two arms. In order to make use of this structure as a 
deterministic entangler, the Bell measurement setup has 
to be replaced through a corresponding projection device 
(post-) selecting that part of the wave function with the 
two electrons distributed between the two arms; alterna- 
tively, this post-selection may be part of the application 
device itself, as is the case in the Bell inequality measure- 
ment. 

In the following, we first derive (Sec. [nj an expres- 
sion for the Bell inequality involving the particle-number 
cross-correlators appropriate for a pulse driven experi- 
ment. We proceed with the calculation of the particle- 
number correlators for a single voltage pulse associated 
with an arbitrary Faraday flux (Sec. IHIfl . The results 
are presented in section Hvl we find the Bell inequalities 
violated for single pulses carrying one Faraday flux, cor- 
responding to one pair of electrons with opposite spin. 
Although the Bell inequality appears to be violated for 
weak pulses (producing less than one pair) too, we argue 
that this violation is unphysical and that its appearance 
is due to a misconception in the original derivation of 



the Bell inequality arising in the weak pumping limit. 
We also generalize the discussion to the situation with 
more complex drives (multi-pulse case and alternating 
pulse sequences) and demonstrate that our Bell inequal- 
ities again are violated only for single-pair pulses flowing 
in either direction through the device. Our analysis of 
an alternating signal produces an apparent violation of 
the Bell inequality, which, however, again appears to be 
an artefact resulting from an improper derivation of the 
Bell inequality for the alternating signal. In both cases 
of failure, weak pulses and alternating pulse sequences, 
we encounter backflow phenomena which spoil the proper 
derivation of the Bell inequality for our setup. 



II. BELL INEQUALITY WITH NUMBER 
CORRELATORS 

The Bell inequality we are going to use here has been 
introduced by Clauser and Horned; it is based on the 
Lemma saying that, given a set of real numbers x, x, y, 
y, X, Y with \x/X\, \x/X\, \y/Y\, and \y/Y\ restricted 
to the interval [0, 1], the inequality \xy — xy + xy + xy\ < 
2|AY~| holds true. We define the operator of electric 
charge Ni(t ac ) transmitted through the i-th spin detec- 
tor during the time interval [0, i ac ], where i ac > is the 
accumulation time. The charge operator Ai(i ac ) can be 
expressed via the electric current Ii{t) flowing through 
the i-th detector, N l {t &c ) = J * ac dt' h{t'). In the Bell 
test experiment, see Fig.^ one measures the number of 
transmitted electrons with a given spin polarization, Aj, 
i = 1, ... ,4, and defines the quantities x = Ai — A3, 
y = N 2 - A 4 , X = Ni + A 3 , and Y = A 2 + A 4 for fixed 
orientations a and b of the polarizers (and similar for x, 
and y for the orientations a and b), see Ref.flrl Our Bell 
setup measures the correlations 



Ky(a,b) = (JV' i (i ac )# J -(t ac )) 



Cat 

■/ 



dhdt 2 (hit^ijih)} 



between the number of transmitted electrons Aj, i 



(1) 



1,3, 



in the lead 'u' with spin polarization along ±a and their 
partners Nj, j — 2,4, in lead 'd' with spin polarization 
along ±b. Using the above definitions for x, y, X, and 
Y , we obtain the normalized particle-number difference 
correlator. 



E(a,h) 



([A!- A 3 ][A 2 - A 4 ]) 
([Ai+A 3 ][A 2 +A 4 ]) 

/C12 — /C14 — /C 32 + /C34 
K-12 + ^14 + JC. 



32 



fC 



34 



(2) 



and evaluating the correlators for the four different com- 
binations of directions a, a and b, b, we arrive at the 
Bell inequality 

E m = \E(a, b) - E(a, b) + E(a, b) + E(a, b)| < 2. (3) 
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We proceed further by separating the current correla- 
tors in Eq. JIJ into irreducible parts Cj,-(a, b; tx, t%) — 
(5ii(ti)5ij(h)) with SIi(t) = Ii{t) - (Ii(t)) and products 
of average currents and rewrite E(a, b) in the form 



£(a,b) 



ifi 



ifi 



K; 



32 



if a 



if 12 + if 14 + if 3 



K- 



34 



A, 



(4) 



where we have defined A± = [{N t ) ± (N 3 )][(N 2 ) ± (AT 4 )] 
with the irreducible particle number correlator 



ifzj(tac) = (^i(tac)^i(tac)) 



dtidt 2 Cij(a.,h;t 1 ,t 2 ) 



(5) 



The average currents are related via (ii(i)) = (h{t)) = 
(7 u (t)>/2and (J 2 (i)) = (J 4 (t)) = (Jd(*))/2 and thus A_ = 
0, A + = (N u )(N d ). The irreducible current-current cor- 
relator factorizes into a product of spin and orbital parts, 
Cy(a, b;ii,t 2 ) = |(aj|b,-)| 2 C ud (*i,t2) with ai >3 = ±a 
and b2. 4 = ±b. The spin projections involve the an- 
gle # a b between the directions a and b of the polarizers, 
(±a| ± b) = cos 2 (6» ab /2) and (±a| =p b) = sin 2 (6» ab /2), 
and the Bell inequality assumes the form 



ifud [cos 6> ab - cos 6> ab + cos 6> ab + cos 6> ab ] 



2A- ud + (A> u )(A> d 



< 1, 



(6) 



where K ud (t ac ) = f Q *° dtxdt^ C u d(*i, h) is the (irre- 
ducible) number cross-correlator between the upper and 
lower leads of the fork. The maximal violation of the 
Bell inequality is attained for the standard orientations 
of the detector polarizations 6* ab = 6* ab = # ab = 7r /4, 
6* ab = 37r/4; the Bell inequality JSJ then reduces to 



2if, 



2K ud + (N u ){N d ) 



< 



1 



v/2' 



(7) 



III. 



NUMBER CORRELATORS FOR A SINGLE 
PULSE 



The orbital part C ud (ti,t 2 ) of the current cross-corre- 
lator between the upper and lower leads can be calculated 
within the standard scattering theory of noissi2ii2i22i2i. 
We assume that the time dependent voltage drop V{t) at 
the splitter can be treated adiabatically (i.e., the voltage 
changes slowly during the electron scattering time). The 
electrons incident from the source lead 's' and scattered 
to the 'up' or 'down' lead then acquire an additional time 
dependent phase <p(t) — f_ dt'eV{t')/%. The scattering 
states (for one spin component) describing the electrons 
in the upper and lower leads, ^ u (x,t) and ~^ d (x,t), take 
the form 



de 



- x/v ^c e + r u a e + t du b e )e 



ikx 




ikx 



-iet/h 



(8) 



b t e~ ikx 



(t^e^-^c, + r d b e + t nA d e )e lkx 



where v e = \/2me; a t , 6 e , and c e denote the annihila- 
tion operators for spinlcss electrons at energy e in leads 
'u', 'd', and 's'; the scattering amplitudes t su (t du ) and 
W (^u d ) describe particle transmission from the source 
(down) lead into the upper lead and from the source (up) 
lead into the lower ('d') lead; r u , r d denote the reflection 
amplitudes into the leads 'u' and 'd'. Such adiabatically 
deformed scattering states (JSJ and © have first been 
used in the calculation of the spectral noise power in an 
ac-driven system 2 ^; the validity of this approach has been 
confirmed in several experiments 2 ^. 

We substitute these expressions into the current op- 
erator I u ( d ){x,t) and drop all terms small in the pa- 
rameter |e — e'|/e F (we assume a linear dispersion). 
The irreducible current cross-correlator C nd (ti,t 2 ) = 
(SI u (x,ti)5I d (y,t 2 )} measured at the positions x and y 
in the leads 'u' and 'd' can be splitted into two terms, 
one due to equilibrium fluctuations, C^(t\ — t 2 ) = 
f{duj/2ir) S^iujy^ 1 -^ with 



S eq {u>) 



2e 2 

■-r- T ud cos(wr+) — 



tijj 



(10) 



and a second term describing the excess correlations at 
finite voltage, 



4e_ 2 
J 2 



^ud(*i)*2j = --To-Ju-'dSin a(7 



with a(r,8) = t: 2 9 2 / smh 2 [k9t /K\ (9 is the temperature 
of electronic reservoirs), r = t\ — t 2 , r ± = (x ± y)/v F , 
£l = t\ — x/v F , and £2 = t 2 — y/v F . The coefficients 
T u = \t su \ 2 , T d = |i sd | 2 , and T ud = |i ud | 2 - |t du | 2 denote 
the transmission probabilities from the source to the 'up', 
'down' leads, and from the 'down' to the 'up' lead. 

The equilibrium part of the current cross-correlator 
C°j(^i — t 2 ) describes the correlations of the electrons 
in the Fermi sea propagating ballistically from lead 
'u' to lead 'd' (or vice versa) with the retardation 
t + = (xi + x 2 )/v F . The corresponding equilibrium 
part of the particle-number cross-correlator, K^i = 

Jq* c dt\dt 2 Cu d (ii — t 2 ) then takes the form 

Kl~ 4 T "dln— , r = max{a/e P) r+} ) (12) 

where we have assumed the zero temperature limit and 
an accumulation time i ac 3> t. The logarithmic diver- 
gence in i ac reduces the violation of the Bell inequality 
Eq. |J7Jl at large accumulation times and one has to sup- 
press the equilibrium correlations between the upper and 
the lower leads in the setup. This can be achieved via a 
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reduction in the transmission probability T u( j, however, 
in the fork geometry of Fig. ^a) the probability T ud can- 
not be made to vanish. Alternatively, one may chose a 
setup with a rcflcctionlcss four-terminal beam splitter as 
sketched in Fig. ^b) with no exchange amplitude be- 
tween the upper and lower outgoing leads; using such a 
fork geometry, the equilibrium fluctuations K™i can be 
made to vanish^. 

Next, we concentrate on the excess part K° d of the 
particle-number cross-correlator (N u (t ac )N d (t ac )) . Note 
that the excess fluctuations are the same for both setups 
Fig. ^a) and (b) and we can carry out all the calcula- 
tions for the fork geometry. We consider a sharp volt- 
age pulse applied at time to, < to < ^ac, with short 
duration St. The total accumulated phase <p{t) then ex- 
hibits a step-like time dependence with the step height 
A<j) = <f>(t + St/2) - <j){t Q - St/2) = -2tt$/$o, where we 
have introduced the Faraday flux $ = — c jV(t)dt and 
$o = he/ e is the flux quantum. The excess part of the 
particle-number cross-correlator K ud then takes the form 
(we consider again the zero temperature limit) 

KS-^J^^^m. (13) 

o 

For a sharp pulse with St <C to,t ac we can identify two 
distinct contributions arising from the integration do- 
mains |*i — *2 1 "C St and \ti — i 2 | ^> St, cf. Refs. l2a 
and [2(1 we denote them with K < and K > . Introducing 
the average and relative time coordinates t = {t\ + t 2 )/2 
and t = t\ — <2 and expanding the phase difference 
<i>{t\) - (f>{t 2 ) = (t>(t + t/2) - <j)(t - r/2) « 0(t)r, the 
first contribution K < reads 



2 ' ac 

= -^r u T d J dt\4>(t)\. (14) 

o 

Assuming that the phase <j)(t) is a monotonic function of 
t (guaranteeing a unique sign for <fr(t)) the last equation 
can be rewritten in terms of the Faraday flux <J>, 

K< = -e 2 T n TA (15) 

<J>0 

this contribution to the particle-number cross-correlator 
describes the correlations arising from the n — 
|$|/<f> additional particles pushed through the fork by 
the voltage pulse V(t), see Eq. l |T8|l below. 

The second contribution K > to originates from 
the time domains < tin) < to — St/2 and to + St/2 < 
*2(i) < t ac , where \<j)(ti) - cj>(t 2 )\ = 2tt$/$ , hence 



here, we have kept the most divergent term in the mea- 
surement time t m = t ac — to, the time during which the 
pulse manifests itself in the detector. The above expres- 
sion describes the response of the electron gas to the 
sudden perturbation V(t); the logarithmic divergence in 
the measurement time t m can be interpreted^ along the 
lines of the orthogonality catastrophe 27 , with the iso- 
lated perturbation in space, the impurity, replaced by 
the sudden perturbation in time. The periodicity of the 
response in the Faraday flux $ is due to the discrete na- 
ture of electron transport as expressed through the bino- 
mial character of the distribution function of transmitted 
particles^SiSi. Remarkably, the above logarithmically di- 
vergent contribution to vanishes for voltage pulses 
carrying an integer number of electrons n — |$|/<!>o, see 
(|18fl below. This follows quite naturally from the invari- 
ance of the scattering amplitudes t su and t s( j in Eqs. ijSJ 
and under the (adiabatic) voltage pulses carrying in- 
teger flux ±n$, t sx — > i sx e ±27rn with x = u, d; trans- 
mitting an integer number of particles at Faraday fluxes 
$ = n$o avoids the system shakeup and the associated 
logarithmic divergence. 

We proceed with the determination of the average 
number of transmitted (spinless) particles (-/V u (d) (/ac)) = 
J* ac dt (I u (d)(x, t)). Within the scattering matrix ap- 
proach the average currents in the upper and lower leads 
are given by the expression 

(4(d) 0>*)) = r T u(d)eF(t- x/v F ) 
c 

= 7^T u{d) (j)(t-x/v F ). (17) 

The time integration provides the average number of 
transmitted particles 

(iV u(d) (t ac )) = er u(d) |-. (18) 

With T u + T d = 1, the result {TBJ tells that a voltage 
pulse corresponding to n — |$|/<f>o flux units pushes n 
spinless electrons through the fork, in forward direction 
from the source lead V to the prongs 'u' and 'd' if <I> > 
and in the backward direction for <f> < 0. 

IV. RESULTS 

Substituting the above expressions for the particle- 
number cross-correlators and for the average number of 
transmitted particles into J7J) we arrive at the following 
general result for the Bell inequality 

= n+(2/7r 2 )sin 2 (7m)ln(W^) (m 
2n 2 -n- (2/7r 2 )sin 2 (7m)ln(W^) ' 1 ' 

A. Pulse with integer flux 



For a voltage pulse with integer n the above expression 
simplifies dramatically as all logarithmic terms vanish, 
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leaving us with the Bell inequality 



Ebi — 



1 



In - 1 



< 



V2' 



(20) 



which we find maximally violated for n = 1 and never 
violated for larger integers n > 1 — any additional par- 
ticle accumulated in the detector spoils the violation of 
the Bell inequality. Furthermore, this violation is inde- 
pendent of the transparencies T u , T<j and hence universal; 
moreover, the Bell inequality (|2U|) does not depend on the 
particular form or duration of the applied voltage pulse 
but involves only the number of electrons n carried by 
the voltage pulse. 

A voltage pulse with n = 1 pushes two electrons with 
opposite spin polarization towards the beam splitter. 
Such a pair appears in a singlet stated and can be de- 
scribed by the wave function 'f^ = 4>l<f>lxlg with the 

spin-singlet state xll = [x\x] ~ x\x}]/V^; (/> s is the or- 
bital part of the wave function describing a particle in the 
source lead 's' and the upper indices 1 and 2 denote the 
particle number. This local spin-singlet pair is scattered 
at the splitter and the wave function transforms to 

where the last term describes two particles in a singlet 
state shared between the upper and lower leads of the 
fork. The Bell inequality test is only sensitive to pairs 
of particles propagating in different arms, implying a 
projection of the scattered wave function ^H at onto the 
spin-entangled component. Thus the origin of the entan- 
glement is found in the post-selection during the cross- 
correlation measurement effectuated in the Bell inequal- 
ity tesljSiS. From an experimental point of view it may 
be difficult to produce voltage pulses driving exactly one 
(spinless) particle n = 1. However, as follows from the 
full expression Eq. 1191) . for a sufficiently small devia- 
tion \Sn\ = \n — 1| <C 1 the logarithmic terms are small 
in the parameter (Sn) 2 and thus can be neglected, pro- 
vided the measurement time t m satisfies the condition 
(cm) 2 ln(£ m /<5£) < 1. 



B. Weak pumping regime 

The weak pumping regime with n < 1, correspond- 
ing to small voltage pulses carrying less then one elec- 
tron per spin channel, deserves special attention. In- 
spection of H19|) shows that the Bell inequality can be 
formally violated in this regime. We believe that this vi- 
olation of the Bell inequality has no physical meaning. 
Below, we will show that for non-integer n many-particle 
effects generate backflow of particles in our setup. We 
argue that this backflow leads to an inconsistency in the 
derivation of the Bell inequality itself, as the assumption 
< \x/X\,\x/X\,\y/Y\,\y/Y\ < 1 may no longer hold 
true for n < 1. 

In order to understand this weak pumping regime bet- 
ter, we analyze the sign of the full current cross-correlator 



C u d(ii,^2) = (Iu(ti)Id(t2)) (for one spin component); ex- 
pressing this quantity through the phase </>(£) we find the 
form 



Cud(^l) £2) — 



(2*) 



<Kti)«K* 2 )-4 



sin 2 ^)-^)/^] 



(ti " £ 2 ) 2 



(21) 



We consider again the case of a narrow voltage pulse 
applied at £ = to and assume a specific shape </>(£) = 
2narctan[(£ — to) /St]. Choosing times £1 < to < £ 2 before 
and after the application of the pulse at to, we find that 
the currents in the leads 'u' and 'd' predominantly flow in 
opposite directions: for a sharp pulse with |£i.2 — to\ 3> St 
we can assume that 4>{t\) — ^(£2) = 27m and thus the 
correlator C ud (£i,£2) takes the form 



C u d(£i, £2) 



7T 2 (5£) ; 



:T u T d 



■ 2 
sin irn 



(1+Z 2 )(1 + Z 2 ) (Z!-Z 2 ) 2 



(22) 



where we have introduced z i:2 — (£1,2 — £0)/ 'St. For 
I #1,2 1 ^ 1 the second (negative) term oc l/(|zi|+z 2 ) 2 de- 
scribing the irreducible correlations dominates over the 
first (positive) term oc l/(z 2 z|) and hence the full cur- 
rent cross-correlator is negative. This negative sign tells 
us that, despite application of a positive voltage pulse 
with n > 0, the currents at times £1 < £0 < £2 in leads 'u' 
and £ d' flow in opposite directions on average. Note that 
this unusual behavior is a specific feature of time depen- 
dent voltage pulses and does not appear for a constant 
dc voltage with </>(£) = eVt/Ti — in this case the current 
cross-correlator is always positive. 

As a consequence, the time-integrated full particle- 
number cross-correlator (per one spin component) may 
turn out negative as well and it does so for voltage pulses 
carrying less then one electron per spin channel n < 1, 



^-ud — e T u Td 



2 • a/ \ t £ 
— — sin (7m) m 



St 



(23) 



Hence, in the weak pumping regime the particles in the 
outgoing leads 'u' and 'd' are preferentially transmitted 
in opposite directions. Note that both (negative) terms 
in the correlator, the one (— n) from short time differ- 
ences as well as the contribution (oc — sin 2 (n7r) ln(£ m /<5£) 
related to the 'orthogonality catastrophe' dominate over 
the (positive) product term (n 2 ), with the second one be- 
coming increasingly important at large measuring times. 
Furthermore, this second term also drives the particle- 
number cross-correlator negative at large non-integer 
n > 1 and long measuring times, again signalling the 
presence of particle backflow in the device. 

The derivation of the Bell inequality relies on the 
assumption that the quantities \x/X\, \y/Y\ etc. are 
bounded by unity. For our setup this implies that the 
particle number ratios of the type \x/X\ = \(N% — 
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N 3 )/(Ni+N 3 )\ are bounded by unity, which is only guar- 
anteed for particle numbers with equal sign, NiN 3 > 0; 
hence, particles detected in the pair of spin filters with 
polarization ±a in the upper arm have to be transmit- 
ted in the same direction. Next, we note that particles 
with opposite spin propagate independently and hence 
our finding that particles preferentially propagate in op- 
posite directions of the outgoing leads 'u' and 'd' for 
n < 1 also implies that the particle numbers N\ and 
A 3 can be of opposite sign, hence the condition x < 1 
is not necessarily satisfied for n < 1. On the contrary, 
for n — 1 the full particle-number cross-correlator (for 
one spin component) vanishes, /C ud = 0: in the simplest 
interpretation we may conclude that the single transmit- 
ted particle is propagating either through the upper or 
the lower lead, thus either 7V U = and iV d — 1 or vice 
versa and the quantity |a3 / -X" j is properly bounded. The 
above arguments cannot exclude the relevance of addi- 
tional many-particle effects, i.e., the appearance of addi- 
tional particle-hole excitations in the system contributing 
to the particle count in the various detectors. A formal 
proof of the unidirectional propagation of particles con- 
firming the applicability of the Bell inequality for the 
present non-stationary situation relies on the calculation 
of the full counting statistics of particles measured in the 
detectors 1 and 3, etc.; such a calculation has not been 
done yet. 



C. Many integer-flux pulses 

Above, we have concentrated on the situation where 
only a single voltage pulse has been applied. Let us 
consider another situation where a sequence of voltage 
pulses driving an integer number of electrons is applied 
to the source lead 's'. In contrast to the previous analy- 
sis, we study the total transmitted charge from t = — oo 
tot — oo, Ni(oo) — dt Ii(t); the excess part of the ir- 
reducible particle- number cross-correlator takes the form 
[we remind that the equilibrium part can be quenched in 
going to the four-terminal beam splitter of Fig. [Jb)] 



Kl = oT u T d / dhdt 



sin 2 [(0fti)-0fe))/2] 
(£i - t 2 ) 2 



(24) 



In our further analysis, we closely follow the technique 
developed in Ref.|23- The double integral in the above ex- 
pression is logarithmically divergent at large times ti,t%, 
producing the logarithmic dependence on the measure- 
ment time t m noted above for the finite accumulation 
time. However, for pulses with an integer number of elec- 
trons, this problematic term disappears; in this case we 
are allowed to regularize the integral in l|24|) with the help 
of 



1 



(*i - t 2 y 



i 



i 



(t x - t a + i8f (h -t 2 - i5f 



(25) 



and 6 — » a small cutoff. Expressing the factor sin 2 (...) 
in Eq. Q24JI in terms of exponential functions, we arrive 
at the form 



A ud — 



e 2 T u T d 
(2tt) 2 



dt\dt2 



e i4>(ti)-i4>(t 2 ) 

(t x - 1 2 + i8f 



(26) 



e i0(ti)— i#(t2) 



In order to proceed further, we split the exponential into 
two terms, e^W = f+(t) + /_(*), with f+{t) and /_(t) 
two bounded analytic functions in the upper and lower 
complex-t plain. Substitution into the above expression 
and using Cauchy's formula for the derivative, 

allows us to write the particle-number correlator in the 
form 

ki = ~^- T ^ J dt [/+(*)/;(*) - /-(*)/*(*)] • (27) 

In (I27|) we have made use of the analytical properties of 
f±(t); in particular, with the complex conjugate func- 
tions (t) and ft (t) bounded and analytic in the lower 
and upper half-planes, respectively, we easily find that 
/ dtf + {t)f*_{t) = J dtf-{t)fX{t) = 0. In addition, we 
can also express the average number of transmitted par- 
ticles in terms of the functions f±(t) introduced above, 



(iV u(d )) = -T n(d) / dtcj>(t) 



(28) 



±T u{A) (dte-^U^ 
iri at 



e 
ni 



rTu(d) / *[/+(*)/;(*) + /-(*)/•(*)]. 



Rewriting Eqs. H27|l and (|28|l in terms of the real numbers 
n±, 



n±=±-L / dtf ± (t)f* ± (t), 



(29) 



we obtain the particle-number cross-correlator and the 
average number of transmitted particles in the form 

K°l = -e 2 T u T d (n++n_), (30) 
(A u(d) ) = 2eT u(d) (n+ - n_). (31) 

Substituting these expressions into the Bell inequality 
Eq. |0, we arrive at the result 



Em 



71+ + n_ 



2(n + — ?i_) 2 — (n_)_ + n_) 



< ^ (32) 



The physical meaning of the numbers n± is easily iden- 
tified for the specific form of Lorentzian voltage pulses 



2^7,/e 



i+it-t^r 



(33) 
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where the index i denotes the number of the pulse in 
the sequence, fcj is the moment of its appearance, is 
the pulse width, and rii the number of spinless electrons 
carried by the i-th pulse with the sign of rij defining the 
sign of the applied voltage. Such a sequence of pulses 
produces the phase 

from which the decomposition into the terms f±(t) can 
be found. The further analysis is straightforward for uni- 
directional pulse sequences with all m > 0, in which case 
exp[i(/>(i)] = f+(t) and n + = J2i n i, n - = 0, or all m < 
whence exp[i(/)(t)] = f-(t) and n + = 0, n- = —J2i n i- 
It then turns oufe^LSS that all results for the irreducible 
particle- number cross-correlator (|30[) . the average cur- 
rents (|31|) . and the Bell inequality (|32|) do neither de- 
pend on the separation tj+i — U between the pulses nor 
on their widths . Furthermore, the result 1)3 2() for the 
Bell inequality agrees with the previous expression 12U|) 
where a single pulse is carrying n — n + (or ) electrons 
in one go and we confirm our finding that the violation of 
the Bell inequality is restricted to pulses containing only 
one pair of electrons with opposite spin. Also, we note 
that for the case of well separated pulses we can restrict 
the accumulation time over the duration of the individ- 
ual pulses, in which case the Bell inequality is violated 
for all pulses with \rii\ = 1. 

Another remark concerns the case of an alternating 
voltage signal with no net charge transport and hence 
zero accumulated particle numbers {N u r^\) = 0. Equa- 
tion 1|31|) then tells us that n + = n_ and the Bell 
inequality (|32[) is formally violated. However, we ar- 
gue that this violation is again unphysical and due to 
the same improper normalization of the basic quantities 
|x/A|, |y/y|, etc. as encountered previously for the case 
of small Faraday flux n < 1: concentrating on the ex- 
pression x/X = (Ni ~ N 3 )/(N 1 + N 3 ), we note that 
two pulses with opposite signs allow for processes where 
the charge driven through the two spin detectors satisfies 
N1N3 < and hence |x/A| > 1, in contradiction with 
the requirements of the lemma. Note that the manner 
of violating the Bell inequality is quite different for the 
physical cases involving pulses with a single particle (see 
the discussion of single integer-flux pulses in Sec. IIV Al 
with n = 1, or the discussion of many integer-flux pulses 
in Sec. IIV CI with n+ = 1, n_ = and n + = 0, n_ = 1) 
and for the unphysical situation of an alternating sig- 
nal with n + = n_ discussed above: in the first case the 
small denominator results from a cancellation between 
the product term 2n and the negative number correlator 
— 1, hence E BI = |l/(2— 1)| = 1, while in the second case, 
the product term vanishes and there is no compensation, 
although the final result is the same, E BI = |1/ — 1| = 1. 
The same apparent violation appears at large non-integer 
values of n > 1 and long measuring times, where the term 
oc sin 2 (7rn) ln(£ m /5f) becomes dominant, cf. I|19|) . 



V. CONCLUSION 

The application of voltage pulses to a mesoscopic 
fork allows to generate spin-entangled pairs of electrons 
through post-selection; the presence of these entangled 
pairs can be observed in a Bell inequality measurement 
based on particle- number cross-correlators. A number 
of items have to be observed in producing these entan- 
gled objects: i) Equilibrium fluctuations competing with 
the pulse signal have to be eliminated. This can be 
achieved with the help of a four-channel beam splitter 
as sketched in Fig. ^b) where the channel mixing is 
tuned such that the transmission T u( j between the up- 
per and lower channel is blocked, it) Pulses V(t) with 
integer Faraday flux <!> = — c j dtV(t) = n$o injecting 
an integer number of particles shall be used. Otherwise, 
the 'fractional injection' of a particle induces a long-time 
perturbation in the system producing a logarithmically 
divergent contribution to the excess number correlator. 
The flux $ = (n + 5n)&o, n = an integer, associated 
with the voltage pulse has to be precise within the limit 
(Sn) 2 -C l/hi(t m /St), with t m the measurement time 
of the pulse and St the pulse width. Hi) The Bell in- 
equality is violated for pulses injecting a single pair of 
electrons with opposite spin, i.e., pulses with one Fara- 
day flux and hence n = 1. The maximal violation of the 
inequality points to the full entanglement of the pair — 
the question what type of pulses produce only partially 
entangled states (as quantified in terms of concurrence or 
negativity of the partially transposed density matrijiSS) 
has not been addressed here, iv) Although weak pump- 
ing with pulses carrying less than one Faraday flux, i.e., 
n < 1, formally violate the Bell inequality (note the pro- 
viso ii), however), we associate this spurious violation 
with an improper normalization of the particle-number 
ratios (Ni — Nj)/(N + Nj) entering the Bell inequal- 
ity, v) The same argument also applies to the case of 
pumping with an alternating signal — we find the Bell 
inequality always violated when the average injected cur- 
rent vanishes (i.e., when the number of carriers transmit- 
ted in the forward and backward directions are equal). 
Again, the origin of this spurious violation is located in 
the improper normalization of the particle-number ratios 
(Ni - Nj)/(Ni + Nj) for this situation. 

The above points suggest the following physical inter- 
pretation: An integer-flux pulse with Faraday flux n$o 
extracts exactly n electron pairs from the reservoir which 
then are tested in the Bell measurement setup. For n = 1 
we find the Bell inequality maximally violated, implying 
that the electrons within the pair are maximally entan- 
gled and not entangled with the remaining electrons in 
the Fermi sea. On the other hand, the application of a 
fractional-flux pulse with non-integer n produces a super- 
position of states with different number of excess electron 
pairs in the fork. The electrons injected into the fork then 
remain entangled with those in the Fermi sea and their 
analysis in the Bell measurement setup makes no sense. 

In our analysis of the spurious violations of Bell in- 
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equalities for weak pumping and for alternating drives 
we have identified the presence of reverse particle flow 
as the problematic element. In the weak pumping limit 
this conclusion has been conjectured from the appearance 
of negative values in the current cross-correlator, imply- 
ing negative values of the particle-number correlator for 
n < 1. Although we believe that these are strong argu- 
ments supporting our interpretation, we are not aware of 
a formal analysis of the backflow appearing in this type 
of systems. The question to be addressed then is: Given 
a bias signal driving particles through the device in the 
forward direction, what are the circumstances and what 
is the probability to find particles moving in the oppo- 
site direction (backflow)? A related problem has been 
addressed by Levito\*2H (see also Ref. l3lTl who has de- 
rived the full counting statistics for the charge transport 
across a quantum point contact in the weak ac pumping 
regime and has identified parameters producing a strictly 
unidirectional flow. The corresponding analysis for our 
system remains to be done. 

A similar scheme for producing spin-entangled pairs of 
electrons has been discussed in Ref. 0, where a constant 
voltage V has been applied to the source lead 's'. In this 
case, the source reservoir injects a regular sequence of 
spin-singlet pairs of electrons separated by the voltage 
time Ty — h/eV; the Bell inequality then is violated at 



short times only. The main novelty of the present pro- 
posal is the generation of well separated spin-entangled 
electron pairs in response to distinguished voltage pulses, 
thus avoiding the short time correlation measurement at 
time scales of order Ty. 

Our mesoscopic fork device produces entangled pairs 
of electrons with a probability of 50 %, i.e., half of the 
single-flux pulses will produce a useful pair with one spin 
propagating in the upper and the other in the lower chan- 
nel. The competing events with both particles moving in 
one channel produce no useful outcome. This is similar 
to the finding of Bcenakker et ali& who derive a concur- 
rence corresponding to the production of one entangled 
pair per two pumping cycles. In how far this represents 
an upper limit in the performance of this type of devices 
or what type of entanglement generators are able to reach 
(at least ideally) 100 % efficiency is an interesting prob- 
lem. 
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